Many-body dynamics of p-wave Feshbach molecule production: a mean-field approach 
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We study the mean-field dynamics of p-wave Feslibacii molecule production in an ultra cold gas of 
Fermi atoms in the same internal state. We derive a separable potential to describe the low-energy 
scattering properties of such atoms, and use this potential to solve the mean-field dynamics during 
a magnetic field sweep. Initially, on the negative scattering length side of a Feshbach resonance 
the gas is described by the BCS theory. We adapt the method by Szymariska et al. [Phys. Rev. 
Lett. 94, 170402 (2005)] to p-wave interacting Fermi gases and model the conversion dynamics of 
the gas into a Bose-Einstein condensate of molecules on the other side of the resonance under the 
influence of a linearly varying magnetic field. We have analyzed the dependence of the molecule 
production efficiency on the density of the gas, temperature, initial value of the magnetic field, and 
magnetic field ramp speed. Our results show that in this approximation molecule production by a 
linear magnetic field sweep is highly dependent on the initial state. 

PACS numbers: 03.75.SS, 05.30.Fk, 34.50-s 



I. INTRODUCTION 

Ultra-cold Fermi gases have attracted considerable 
interest due to the possibility of studying many-body 
fermionic physics in well controlled and tunable environ- 
ments Some of this interest has been stimulated by 
the prospect of observing pairing phenomena similar to 
those observed in superconductors and liquid "^Hc. The 
existence of magnetically tunable Feshbach resonances 
in ultra cold gases allows the interaction strength 

between the atoms to be tuned such that the system be- 
haves as a gas of long range Cooper pairs described by 
the theory of superconductivity presented by Bardeen, 
Cooper and Schrieffer (BCS) Q in the limit of weak in- 
teratomic interaction, and as a Bose-Einstein condensate 
(BEC) of diatomic molecules in the limit of weak repul- 
sion Q. The transition between these limits is often re- 
ferred to as the BCS-BEC crossover 

Identical fermions in the same internal state are for- 
bidden from colliding in the s-wave (^=0, £ being the 
value of the relative angular momentum quantum num- 
ber) so that the lowest partial wave scattering amplitude 
becomes the p-wave (^=1). Pairing between atoms in- 
teracting in the p-wave has been theoretically studied 
in the context of '^He and ultra cold gases [Tol-fl^ . 
p-wave Feshbach resonances have been observed in ul- 
tra cold gases of the fermionic isotopes '^^K [l^ and ^Li 
[13, In the case of experiments have shown the 
existence of a splitting of the p-wave resonance which 
was explained by the magnetic dipole -dip olc interaction 
between the atomic valence electrons [1^ . 

The formation of weakly bound diatomic p-wave 
molecules has been experimentally observed in both ''^K 
and "^Li [l3,[il]- In the work on ''"K, atoms were pre- 
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pared in the state labelled |/, m/) = |9/2,— 7/2), where 
/ denotes the hyperfine state of an atom in the absence 
of a magnetic field and m/ denotes the Zeeman state 
of the atom. Molecules were formed using a resonantly 
modulated magnetic field as first demonstrated for 
the formation of molecules in Bose gases [l^. Binding 
energies of the ^''K2 molecules were measured on both 
the BEC side where the atoms form a bound state and 
on the BCS side where the atoms form a resonance state 
localised by the centrifugal barrier. In the latter case 
lifetimes for decay from the resonance state were also 
measured. Similar experiments in ^Li used linear sweeps 
of the magnetic field to form molecules [l^ [2^ [2l| . 

Motivated by these experiments we study the dynam- 
ics of molecule formation from a gas of fermionic atoms 
in the same internal state, which occurs during a linear 
sweep of the magnetic field. Such a system has been stud- 
ied previously by Szymariska ct al. j22| , who used a mean- 
field approach to investigate s-wave pairing for fermions 
in different spin states. Here we extend this treatment to 
deal with the more complicated p-wave pairing case. The 
initial state is taken to be a BCS paired gas on the neg- 
ative scattering length side of the resonance. We study 
the dependence of the molecule production efficiency on 
the speed of the magnetic field sweep, as well as on tem- 
perature, density and initial magnetic field position. We 
also study the evolution of the order parameter and the 
molecule density after a rapid variation of the magnetic 
field. 

This article is structured as follows: In Sec. II, we 
introduce the pseudo-potential that we use to model the 
two body physics in the many body problem. In Sec. 
Ill, we present our BCS model for studying the initial 
state of the gas from which to study the dynamics, while 
Sec. IV contains the dynamical equations. We study 
linear sweeps of the magnetic field from the BCS side of 
the resonance into the BEC side where we calculate the 
molecule production efficiency by taking the overlap of 
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the pair function with the two body bound state. We also 
analyze the prospects to observe atom-molecule coherent 
oscillations as a function of the time after a fast sweep 
through the resonance. Finally, Sec. V is a discussion of 
our results. 



II. TWO-BODY MODEL 

In order to study the many-body effects present in the 
formation of p-wave Feshbach molecules as a linear mag- 
netic field sweep is applied across the resonance, we ex- 
tend the methods of Szymahska et al. [1^ to the case 
of a p-wave resonance. We must therefore choose a suit- 
able representation for the two-body physics. Due to 
the low energies involved in ultra-cold collisions, it is not 
necessary to use the exact form of the interatomic po- 
tential, provided that all relevant low energy scattering 
properties are reproduced. Previous work on fermions 
has been successful in describing the two-body physics of 
Feshbach resonances using pseudo-potentials in both the 
s-wave and p-wave [l^ . Both the mean- field thermody- 
namics and dynamics have been studied using a separable 
potential to model the two-body interaction close to an 
s-wave Feshbach resonance [13, H^l- Motivated by this 
we seek a similar separable form of the p-wave potential 
that can be used in the many-body calculations. 

The presence of a magnetic field defines a quantisa- 
tion axis. Particles interacting in the p-wave have one 
quantum of relative angular momentum and there are 
three possible projections of this vector onto the quan- 
tisation axis. The situation can be encapsulated in a 
three-component separable potential 

V = |xi.i)Ci.i(xi.i| + |xi.o)6,o(xi.o| + |xi -i)6,-i(xi 

(1) 

Here \xe.me) is the form factor of the potential and ^i^mi 
is the amplitude, where denotes the projection of this 
angular momentum onto the magnetic field, which we 
choose to be in the z-direction. This pseudo-potential 
can account for the observed splitting of the resonance 
feature [l6j into distinct resonances depending upon the 
value of |mf|, due to magnetic dipole-dipole interactions 
between the valence electrons. For each individual term 
a suitable choice must be made for the form factor to 
recover the low energy scattering behaviour of the atom 
pair. A model that reproduces the threshold behaviour 
should be sufficient for our calculations since this is the 
energy range in which the experiments have been per- 
formed. 

We set the parameters of the pseudo-potential to repro- 
duce physical observables. One feature of the two body 
physics that we can compare with experiment is the mea- 
sured binding energy of the molecule. Bound states are 
associated with a pole in the T-matrix given by the 
Lippmann-Schwinger equation 



Here Go{z) = {z — Ho)~^ is the free Green's function and 
Hq is the free Hamiltonian. Experimental measurements 
show that the resonances due to the different |m^| values 
in ^'^K are narrow and well separated [l6j . and so we 
assume each term of Eq. ([T]) can be treated individually 
leading to a T-matrix for each component in the form 



|Xlmi)Clmi(x 
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(3) 



A pole in Timi (z) coincides with a vanishing de- 
nominator in Eq. ([3]). The T-matrix is related 
to the scattering amplitude through fi,mi{p) — 

— 7rm?i(p£, jT^^TOj, 4- iO^ |p£, m^), where p is the 

magnitude of the relative momentum and m is the single 
particle mass. Given that the long range behaviour of 
the interatomic interaction is dominated by the van der 
Waals potential [i^, the parameter S,imi can be related 
to the scattering properties of the system from the low- 
energy limit of the partial wave scattering amplitude for 
the s-wave and p-wave j26l | 
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(4) 



Here ai.m^ is the partial wave scattering "length" , which 
in the case of the p-wave has the dimensions of volume 
and shall here be referred to as the scattering volume. 
Furthermore, in the vicinity of a resonance the scattering 
volume can be parameterised by the magnetic field using 
the well known resonance formula [23, Hi 
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Here a^^^ is the partial-wave background scattering 
"length" for the m^th component, ABe^mt is referred to 
as the width of the resonance and B^ is the resonance 
position. 

In Appendix [X] we discuss how we obtain a suitable 
p-wave form factor, and give values for the parameters of 
our separable potential for the case of ^"K atoms. Using 
this potential, we can solve the equation for the bound 

state 



1 - Clmi(Xlmi|G'o (£■-!) IXlmi) = 0, 



(6) 



where is the energy of the least bound state of the 
atom pair. We then get the low energy expansion of the 
bound state energy 



mai 



(7) 



T{z)^V + VGoiz)T{z). 



(2) 



where aimi is a parameter associated with the range 
of the potential, as explained in Appendix El This ex- 
pression shows that the bound state energy varies as 
the inverse of the scattering volume, and to this ex- 
tent is in agreement with previous two-channel calcula- 
tions [28l - l30| . A plot of the bound state energy about 



3 




B(G) 

FIG. 1: (Color online) Variation of the p-wave binding energy 
with the magnetic field for the mi = resonance in '^'^K. The 
solid blue line is the solution to Eq. HA3|) . The dashed red line 
is the low energy expansion of Eq. The black solid line 
shows the variation of the scattering volume (right axis) in 
the vicinity of the resonance. 
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FIG. 2: (Color online) p-wave elastic scattering cross section 
for ^"K colliding in the |9/2, —7/2) channel as a function of 
collision energy at three different magnetic field strengths. 
The solid blue line corresponds to the pseudo potential model. 
The red crosses are from a coupled-channel calculation. 



the Bi g =198.85 G toi = resonance in '^^K found from 
Eq. (|A3|) is given in Fig. [TJ Also plotted is the low en- 
ergy expansion of Eq. 1^ , vifhich shovifs a good agreement 
close to the resonance. 

For p-wave resonances, it has been observed that the 
magnetic dipole-dipole interaction plays an important 
role, leading for example to the splitting of the mi — 
and nil = il resonances [l^. For such an interac- 
tion, which decays as 1/r^ at long distances, no scat- 
tering length can be properly defined However, 
a scattering volume model may give the correct scat- 
tering cross section behaviour within some finite range 
of collision energies. Physically this can be regarded 
as truncating the interaction at the radius where it be- 
comes small compared to the kinetic energy. To illustrate 
this Fig. [2] shows the scattering cross section given by 
our pseudo-potential, and compares it to a much more 
detailed coupled-channels calculation. Here we obtain 
the coupled-channels results using the Manolopoulos log- 
derivative propagation technique (3l| including a dipole- 
dipole interaction term in the Hamiltonian [32|; we use 
the Born-Oppenheimer potentials of ref. [11] . It can be 
seen that our separable potential reproduces the low en- 
ergy scattering behaviour to a good degree over several 
orders of magnitude, particularly close to the resonance, 
and we therefore proceed in the following sections to use 
this potential in the many-body problem of molecule for- 
mation. 



superconductivity Q. While trapped gases usually have 
anisotropic and inhomogeneous density distributions, for 
simplicity we consider here a spatially homogenous sys- 
tem. 

In the following we review some standard results of 
BCS theory |3J| to present the notation we use in the 
following paragraphs. The many body Hamiltonian for a 
system of fermions in the pairing approximation is given 
by 



H 



Y.m\j)ala, 



(8) 



+ i ^ (kllVlmn) (^(aj,a]')a„a„ -|- aj^a| (a„a„ 

klmn 

Here the and a| are the usual fermion annihilation and 
creation operators, which obey the usual fermion anti- 
commutation rules, K is the single particle kinetic energy 
operator, and V is the interaction potential. The brack- 
ets (...) represent an average over the thermodynamic 
state of the system in the grand canonical ensemble where 
particle number is not fixed, \i) are single particle states 
and are two particle states. The finite temperature 
Green's functions of the system can be defined as 



grsir, 



-(T.K(r)at(r')]), 



(9) 



III. MEAN-FIELD THERMODYNAMICS 

The stationary properties of a weakly attractive ultra- 
cold Fermi gas can be described by the BCS theory of 



FUr,T') = -{TAaUT)aliT')]), 



(10) 



the latter representing pairing in the gas. Here TV is the 
imaginary time orde ring operator that puts the smallest 
value of r to right [SJ]. It is also often convenient to 



4 



define a gap function 



A, 



ij 



[a^aj){rs\V\ij), 



(11) 



which acts as the order parameter of our system. 

By working in the momentum representation and con- 
sidering a translationally invariant system we can write 
down the Hciscnbcrg equations of motion for the Green's 
functions. The resulting equations can be Fourier trans- 
formed and combined with Eq. (fTTj) to give 



A*(p) 



A(p|^|q)4r^tanh(/?^), (12) 



constant. Equations p2l) and ((131) ^^'^ often referred to 
as the gap equation and the density equation, respec- 
tively. These are the BCS equations for the system and 
must be solved simultaneously. Experiments are usually 
performed at constant temperature, which fixes /? in our 
system. This leaves the chemical potential, /i, and A(j3), 
as parameters to be solved for. 

Using the potential of Eq. ([T]) the order parameter can 
be written as 



(14) 



and 



V 2 ' 



i-At.„M3^) 

,1/2 



(13) 



where eq = (S^ ^ |A(q)|2)v^^ ^ £. _ V is the 
volume of the system, /? = jp-rp and ks is Boltzmann's 



where A^^ and XmAp) = (pl: "^<?lxi,mj) represent com- 
ponents in the p-wave. We now drop the label £ in the 
subindex in the functions we have defined, since we will 
only be concerned with the p-wave, i = 1. By inserting 
this into Eq. p2|) we arrive at a set of coupled equations 
for all components of the gap in the angular momentum 
basis 



A* 



Xmi {q)Yu,n, (q)Cmi a;;, Xrn[ {q)Y*„^,^ (q) 



2ec 



tanh 



(15) 



Experiments have shown the toi = 1 and mi = — 1 states 
to be degenerate [l^. The degeneracy of the two com- 
ponents is due to the rotational symmetry of the system 
about the magnetic field axis. The parameter A^^ for 
the mi = 1 and toi = — 1 components can be written as 



A 



±1 



±iA. 
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(16) 



where the magnitudes of these components are identical. 

It is possible to reduce the set of three coupled nonlin- 
ear equations to two [35| . so that we can write 



a: 



('0x(q),'0x(q)) (V'x(q), V'^(q)) 
(V'z(q),-0a;(q)) (V'z(q), V'z(q) 



Here brackets (•, •) represent scalar products of two func- 
tions and should not be confused with a thermodynamic 
average. The functions are given by. 



V'x(q)=n/^^^^^Xii(9)sin0cos^, (18) 



An 



and 



V'z(q) = n/ ^A^^xo(g) cos 6*, 



(19) 



where for brevity we have defined 



/(q) 



tanh 



2^q 



(20) 



We have solved Eq. (fTT]) for the resonance around 
198.85 G in -^^K. We found that the off-diagonal ma- 
trix elements are several orders of magnitude smaller 
(17) than the diagonal elements for densities of 10^"^ cm~^ 
to 10^^ cm"'^. In this case neglecting the off-diagonal el- 
ements does not noticeably change the solution. A plot 
of the solutions to the BCS equations for "'"K is shown 
in Fig. [3l for the gap parameters, A^^ and Fig. |4] for 
the chemical potentials fi. Previous work has used the 
Bose-Fermi 36 , 37. ] model to study the thermodynamic 
properties of resonant p-wave gases p^l - [l2i \3§L l39j . We 
have calculated the chemical potential and gap parame- 
ters using both models and shown the results of the two 
models to be very similar throughout the crossover re- 
gion [35j . 
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FIG. 3: (Color online) Variation of the parameter Am with 
magnetic field for the p-wave resonance in ''"K for a density 
of lO^'^cm"^ and a temperature of 70nK. The solid green line 
is the value of the m = resonance and the dashed blue line 
is for the |m| = 1 resonance. There is no significant differ- 
ence between the results obtained when coupling between the 
components is included and when the coupling is excluded. 
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FIG. 4: (Color online) Values of the p- wave chemical potential 
for the resonance in ^''K for a density of lO^^cm"^ and a 
temperature of 70nK. The dashed blue line is the |mi| — 1 
resonance and the solid green line is the mi — resonance. 
The values for coupled resonances and separated resonances 
are indistinguishable. The positions of the resonances are 
marked by the vertical lines. 



IV. MEAN-FIELD DYNAMICS 

A. Dynamics of molecule production 

Experiments on the formation of molecules often use 
finite speed magnetic field variations to create molecules 



from gases of ultra-cold fermions 0, |40| - |4^ . This has 
stimulated theoretical research into studying the dynam- 
ics of molecule production in these gases [22, 143l - l50l | . At 
the mean-field level the gas can be described by two corre- 
lation functions: the pair function, = {ajai)t, and 
the one body density matrix Tij{t) ~ (a|ai)t. The brack- 
ets (...)t represent a thermal average at time t. Their 
Heisenberg equations of motion in the momentum repre- 
sentation for a homogeneous system are [l^l 

ih^r{p,t) = 2(2^;^)3/2^ Im($*(p,t)(p|V|$(t))) (21) 

ih-^^^p,t) = {p\H2Bmt)) 

-(p|y|$(i))r(-p,i)(2^7i)'/2 

-(p|y|$(i))r(p,t)(27r?i)3/2. (22) 

where due to translational invariance we have Fp^p = 
r(p)(5(0) and $p,-p = <I'(p)(5(0). The initial correlation 
functions are related to the solutions of the BCS equa- 
tions of Sect, mil through 

and 

A(g) - J d\ (p|y|q)ci>(q) (24) 

These equations form a closed set of equations for the 
density distributions, without need for auxilary param- 
eters or constraints p^ . Dynamical equations can be 
written for the partial wave components of both the pair 
function and the one-body density matrix. For example, 
the component of the pair function is given by 

^i,mM^t) ^i' j dn y,;„^ (r!)$(q,t), (25) 

where 17 is the solid angle in q. A similar expression can 
be written for the density matrix components, F^.m^ (g, t). 
This allows the angular integrals in the dynamic equa- 
tions to be performed analytically, reducing the dimen- 
sionality of the numerical problem provided the contribu- 
tion from higher partial waves is not significant. When 
the initial conditions allow the resonances to be treated 
separately the initial state of the gas is in either an 
|mi| = 1 state or an mi = state. At this level of ap- 
proximation it is not possible to transfer population be- 
tween these two components through the dynamics and 
we therefore treat them as individual cases. We study 
the case of a linear magnetic field sweep 

B{t) = Bi- Bt (26) 

where Bi is the initial magnetic field above the resonance, 
i = is the time at which the initial state is prepared 
and B is the sweep rate. At some point the sweep will 
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FIG. 5: (Color online) Variation of the molecular production 
efficiency in as the initial value of the magnetic field is var- 
ied. The different lines represent different densities. Higher 
densities give a larger production efficiency. The temperature 
is 70nK and the sweep rate is held constant at B = 60 G/ms. 
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FIG. 6: (Color online) Variation of the molecular production 
in *''K efficiency as the atomic density is varied. The sweep 
rate is held constant at i? = 60 G/ms. The lines represent 
temperatures of 70nK (solid blue) , lOOnK (dot-dashed green) 
and 200nK (dashed red). 



cease to produce molecules, so that we may choose a final 
magnetic field in the BEC side of the resonance where the 
molecule production has saturated. 

The density of molecules at the end of the sweep can 
be calculated from a Wick expansion of the two-body 
correlation function. This can be approximately written 
as an overlap with the two-body bound state [22| 



1 



d^p 0:(p)$(p) 



(27) 



where (p) is the two-body bound state wavefunction 
in the momentum representation obtained by solving the 
Schrodinger equation using the separable potential. In 
Eq. (P7)) the pair function, $(p) is evaluated at the final 
magnetic field strength of the sweep. 

As a first approximation we retain the lowest order 
partial waves. The effect of adding in the higher partial 
waves can then be investigated. To study the dynamics 
of mi =0 molecule formation in this approximation we 
retain the functions $io(p, i) and roo(p, i). 

We have investigated the variation in mi = molecule 
production as a function of the initial magnetic field, 
In this case the temperature, atomic density and ramp 
speed of the magnetic field are held constant. The re- 
sult is plotted in Fig. [51 As the initial magnetic field 
moves closer to the the resonance position the molecule 
production increases. This is due to an increase in the 
pairing in the initial state of the gas, as can be seen 
from Fig. [31 Similarly, increasing the atomic density in- 
creases the production efficiency as shown in Fig. [6l We 
remark that our mean-field calculations predict a satu- 
ration of the molecule production efficiency at increasing 
density. This is an expected many-body effect, in con- 



trast to purely two-body calculations that would predict 
Nm/Na to be linearly proportional to the density. 

Investigating the variation of the production efficiency 
as a function of initial magnetic field proved to be difficult 
at low densities (^ 10"'^^ — 10^'^ cm"'^). This is because 
the gap parameter quickly drops to zero on the BCS side 
of the resonance at these low densities. At such densities 
it should then be possible to approximate the production 
efficiency at the mean-field level with an instant projec- 
tion of the initial pair state onto the bound state. This is 
because when the initial field is close to the resonance the 
quantities of interest will not evolve significantly before 
the final magnetic field is reached. 

It has been shown by Iskin et al. [5l| that for a har- 
monically trapped p-wave superfluid the central density 
is much larger than for the corresponding s-wave super- 
fluid. This indicates that it may be feasible to reach 
higher densities in p-wave Fermi gases. Following this 
observation, we performed calculations setting the initial 
atomic density to a value well above those normally used 
in experiments (greater than 10^^ cm""^). In Fig. [7l we 
present calculations performed for a initial atomic den- 
sity of 10^^ cm~'^. These results show that it is possi- 
ble to produce a significant number of molecules even 
when pairing in the initial state of the gas is not very 
large. This can be achieved by adjusting the sweep rate 
to lower and lower values, thus allowing for two-body 
correlations to build up during the sweep via intrinsi- 
cally many-body processes. The closer the value of Bi is 
to the resonance the more difficult it is to influence the 
number of molecules produced by varying the sweep rate. 
This is because there is already a great amount of pairing 
in the gas and the production efficiency Nm/Na can not 
exceed 0.5. On the other hand, at low ramp speeds the 
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FIG. 7: (Color online) Variation of the molecular production 
efficiency in as the ramp speed is varied for a temperature 
of 70nK. The blue, dot-dashed curve refers to an atomic den- 
sity of 10^^ cm~^ and an initial magnetic field of 199 G. The 
green, dashed curve refers to an atomic density of 10^^ cm~^ 
and an initial magnetic field of 200 G. The red, solid curve 
refers to an atomic density of 4 x 10^* cm~^ and an initial 
magnetic field of 200 G. 



production efficiencies all converge and become indepen- 
dent of following the adiabatic value obtained from 
the thermodynamic state. 

It may still be possible to produce molecules from a 
linear sweep of the magnetic field at lower densities. As 
a comparison, the molecule production at atomic den- 
sities of 4 X lO^^cm^"^ as a function of inverse sweep 
rate has been plotted in Fig. [7] (solid, red line). The 
initial magnetic field has been set to 200G so that a 
direct comparison can be made with the dashed-green 
line. More molecules are produced relative to the produc- 
tion from an infinitely fast sweep for the lower density of 
4 X lO^^'^cm"'^. However, to achieve a similar net produc- 
tion efficiency as the lO^^cm"'^ gas, the ramp speed will 
have to be set to significantly lower values where eff'ects 
due to higher order correlation functions (not captured in 
the present mean-field approach) may become relevant. 

To analyze how many molecules are actually produced 
during the dynamics, we show in Fig. [8] the fraction of 
molecules produced from a 500G/ms sweep subtracted 
from the fraction produced from a lOG/ms sweep. It 
can be seen that there is a maximum in the production 
efficiency at a certain value of the density. This peak in- 
creases in magnitude and moves to higher densities as the 
initial value of the magnetic field Bi is moved to higher 
fields. This shows that combinations of initial density 
and initial magnetic field can be picked to optimise the 
production efficiency. 

We have checked that including higher partial wave 
components to study the dynamics does not significantly 
affect the production efficiencies reported here |35| . It ap- 
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FIG. 8: The difference in the molecule production efficiency 
from a sweep of lOG/ms and a sweep of 500G/ms as a function 
of density. Here, is the number of molecules over the 
number of atoms after a sweep at a speed equal to B in G/ms 
for a temperature of 70nK. This shows how many molecules 
are actually produced during the dynamics. The dashed curve 
refers to an initial magnetic field position of Bi — 199. 3G, 
while the lower curve refers to Bi — 199. 2G. It can be seen 
there is an optimum density at which to produce molecules 
from the dynamics. 

pears that retaining only the lowest £ components is suffi- 
cient to calculate the production at the mean-field level. 
Similar results are obtained for the mi = 1 molecules 
by allowing the initial state of the gas to be paired with 
TOi = 1 symmetry. 



B. Atom-molecule coherence 

In experiments on s-wavc molecule production, rapid 
sweeps of the magnetic field were used to probe the state 
of the Fermi gas in the region about the resonance. It was 
hypothesised that if the magnetic field was swept into the 
BEG side fast enough, such that the typical sweep time 
was less than the typical collision time, then it would 
be possible to extract information about the gas in the 
strongly interacting region [s^ . The question then arises 
of how the state evolves after such a sweep. If the final 
state, held at a fixed field value, undergoes processes that 
significantly change it, then this method may not be a 
reliable way of probing the gas. For the s-wave it has 
been shown that under such a magnetic field variation the 
final molecule production efficiency will oscillate but with 
a small, decreasing amplitude [l^ . We use an essentially 
identical method to show that this is also true in the 
p-wave and it would not be possible to observe atom- 
molecule coherence with this approach. 

Fig. [5] shows the variation in the production efficiency 
as a function of time after such a magnetic field variation. 
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In this figure, the difi^erent fines correspond to difi'erent 
final magnetic fields. The variation in the molecule pro- 
duction over this time period is given as a percentage and 
seen to be on the order of 10~^ %, which is very small. 
The oscillations in the production are heavily damped 
with the oscillation frequency and damping rate increas- 
ing as the final magnetic field moves away from the res- 
onance. For the case where the final field is located at 
196.5 G the oscillations are not visible on the scale of the 
figure after 20 fis. 

In Fig. [To] the initial magnetic field is varied and the 
final magnetic field held constant. Again the oscillations 
for all detunings are on the order of 10""^%. Both the fre- 
quency and amplitude of the oscillations increase as the 
initial field moves further from the resonance, but not sig- 
nificantly. It should be noted that this appears to be in 
contrast to the s-wavc where the amplitude increases as 
the initial field moves towards the resonance [22|. How- 
ever, in the s-wave studies the final field was generally 
chosen to be much further from the resonance position, 
and so corresponds to a BCS like state projected quickly 
onto a deeply bound BEG state. The near orthogonality 
between the initial and final states in such a case means 
that oscillations are likely to be small in the mean-field 
limit. In our results with the final field closer to the res- 
onance position we see merely the change in the pairing 
between initial and final state reflected. In both cases 
the amplitude of the oscillations is very small (the re- 
sults compared to in the s-wave correspond to a density 
of 1.5 X 10" cm-3). 

We have also studied how the order parameter varies 
after such a magnetic field variation. In this case the gap 
parameter is a function of time defined by 



A(i) = ? / A(x|q>(q|$(i)>, 



(28) 



where wc have used the separable potential to divide out 
a form factor from each side of the equation. We note 
that the value of the binding energy does not enter this 
equation directly. We compare this value against the 
value of the gap parameter when the system is in equi- 
librium at the final magnetic field position. We note that, 
in general, the quantity in Eq. p8)) is complex. As for the 
case of the density variation, we have studied the effect 
of varying both the initial and final magnetic fields. 

We plot the time evolution of the gap parameter in 
Fig. [TlJ The top and bottom panels refer to final mag- 
netic fields of 197G and 198G, respectively. In each panel 
the different lines correspond to different initial mag- 
netic field positions of 198. 2G (top, red), 200. 2G (middle, 
green) and 201. 2G (bottom, blue). It can be seen that 
the closer the initial and final fields are to each other the 
closer the value of the gap parameter is to the stationary 
state value at the final magnetic field position, denoted 
here by A^q. In all cases, the oscillations have a small 
amplitude and quickly decay. The real and imaginary 
parts of the gap parameter oscillate at a frequency that 
corresponds to the energy of the bound state at the final 
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FIG. 9: (Color online) Evolution of the molecule production 
efficiency after an infinitely fast sweep of the magnetic field 
across the 198.85 G resonance in '"^K. The initial magnetic 
field is 199 G, just above the resonance. The different lines 
correspond to differing final magnetic fields of 198.5 G (sofid, 
green line), 197.5 G (dashed, blue line) and 196.5 G (dot- 
dashed, red line). n{t) is the density of molecules as a function 
of time where n(0) is the density of molecules directly after 
the magnetic field variation. 



magnetic field position. This is expected and serves as a 
test on the numerics. The value of oscillation frequency 
of the absolute value of the gap parameter increases as 
the initial magnetic field moves away from the resonance 
position and approximately corresponds to the sum of 
the final bound state energy and twice the initial chem- 
ical potential energy, a measure of the change in energy 
of the paired state during the sweep. 

We conclude that it would not be feasible to observe 
atom-molecule oscillations in this p-wave resonance due 
to the small, rapidly decaying amplitude of the den- 
sity oscillations. This is essentially the same conclu- 
sion reached in Szymanska et al. |22| but extends this 
result to the closed-channel dominated p-wave resonance 



in "K. This also suggests that the method of fast sweeps 
to probe a fermionic p-wave paired condensate would be 
a suitable method to probe the condensate were such 
conditions favourable. 



V. CONCLUSIONS 

We have studied the dynamics of p-wave Feshbach 
molecule formation by linear sweeps of the magnetic field 
within a mean-field approach. It turns out that the 
molecule production efficiency is highly dependent on the 
initial state of the gas. As the density increases it is in- 
creasingly possible to explore a range of initial magnetic 
field values, Bi. As the value of Bi moves away from the 
resonance and deeper into the BCS region it becomes 
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FIG. 10: (Color online) Evolution of the molecule production 
efficiency after an infinitely fast sweep of the magnetic field 
across the 198.85 G resonance in *°K. The final field is held 
constant at 197.5 G. The different lines represent different 
initial magnetic fields of 199 G (solid, blue line), 199.5 G 
(dashed, green line), 200 G (dot-dashed, red line) and 200.5 G 
(dotted, black line). Other symbols are as defined in Fig. [9] 



easier to produce molecules through a sweep of the mag- 
netic field if the density is high enough. Moreover, for a 
given magnetic field value there exists a density at which 
molecule production from a linear sweep is optimal. This 
is due to the fact that as the density increases so does the 
value of the molecule production from an instantaneous 
projection of the initial state onto the final molecular 
wave function. We have treated the mi = and |mi| = 1 
molecules separately since even at relatively high density 
the thermodynamic states are well described separately. 
In turn this allows the dynamics to be treated as if the 
resonances were uncoupled. The initial assumption that 
the partial wave series for the dynamical functions will 
converge seems to be justified since adding higher partial 
waves has little effect on the production efficiency. We 
have also shown that, similarly to the case of an s-wave 
paired system it would be very difficult to observe 
atom-molecule coherence as the result of a fast sweep 
experiment. 

We note that our production efficiencies are generally 
reasonably high, and in our treatment a sizeable fraction 
of the atoms can be associated into molecules using ex- 
perimentally realistic magnetic ramp speeds. However, 
we have not included relaxation processes into our treat- 
ment as such terms are beyond the mean-field theory we 
have used. Experiments to date on p-wavc molecules have 
found them to be short lived, with lifetimes of between 
2 and 20 ms [H, [13, [llll^l , due to dipolar relaxation 
processes [l3l and three-body losses [H, [H, [IJ. Our 
results indicate that small experimental molecular effi- 
ciencies cannot be explained by mean-field theory. 
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FIG. 11: (Color online) Variation of the quantity |A(f)|/Aeq 
with time for final magnetic fields of 197 G (top panel) and 
198 G (bottom panel). The different lines correspond to dif- 
ferent initial magnetic field positions of 198.2 G (top, red), 
200.2 G (middle, green) and 201.2 G (bottom, blue). 
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Appendix A: Modelling of the two-body interaction 

Previous work has successfully used separable poten- 
tials with form factors of a Gaussian form to model the 
s-wave two-body interaction in the context of Feshbach 
resonances [12, . This is convenient both analytically 
and numerically because Gaussian integrals can often be 
performed analytically. We wish to describe the system 
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\me\ 


B,V^(G) AB(G) 







198.85 -21.9482 


-1049850 


1 


198.373 -24.9922 


-905505 



TABLE I: Parameters used for modelling the two-body inter- 
action, (as indicates the Bohr radius.) 



close to threshold and therefore choose the form factor 
to reflect the physics in this region. In the p-wave we 
choose the form factor to be 



(PiXr, 



■kTi 



b/2 



"yi,"H(P) 



(Al) 



which ensures the correct threshold behaviour of the 
bound and scattering states Here is a length 

scale associated with the range of the potential and 
Yi^mi(p) is the £ = 1 spherical harmonic with a projec- 
tion TOi onto the chosen z-axis, in this case the magnetic 
field axis. 

Using the form factor above and the low energy limit 
of the scattering amplitude it is possible to relate the 
scattering volume to the parameters of the separable po- 
tential amplitude, ^„ij defined in Eq. ([1]), and the form 
factor, Eq. 



air. 



(A2) 



where x = m^mi /(47r?i^(Tmj ) is dimensionless, with m 
being the mass of the atom. Equation ^ for the bound 
state energy is then given by 



1 



1 



l-\/¥yey'erfc(y)]) =0, (A3) 



(< 0) 
exp(- 



is the bound 

u^) du is the 



where y = ami\/-^rnE~i/h, E- 
state energy, and erfc(2/) = ^ 
complementary error function. 

All of our calculations have been performed for the 
fermionie isotope ^°K. We have fixed the parameters of 
our separable potential using the experiiTiental results of 
Refs. |l6l [13]. To do this, we take the low energy ex- 
pansion of the binding energy as given by Eq. ([7]) and 
parameterise the scattering volume as in Eq. ([5|). Close 
to the resonance we can expand the resonance formula as 
a series in powers of {B — the terms of which can 

be directly read off from the expression given in Ref. [l^ 
fixing the resonance parameters in Eq. ([S]). The values 
obtained are given in Table IH 

Equation ([7]) can be differentiated with respect to 
the magnetic field to give the magnetic moment of the 
molecules relative to the free atoms. This can be di- 
rectly equated to the experimental values of the magnetic 
moment allowing the parameter ami to be calculated. 
Ref. [13 gives the magnetic moments of the molecules as 



and 



dE 
'dB 

dE 
'dB 



m-i— 



188±2 kHz/G, 



193± 2kHz/G. 



(A4) 



(A5) 



|mi| = l 

The corresponding range parameter is given by 



m.ABr, 



2 "'me 



dE 
'dB 



(A6) 



[1] S. Giorgini, L. P. Pitaevskii, and S. Stringari, Reviews of 

Modern Physics 80, 1215 (2008). [14] 
U. Fano, Nuovo Cimento 12, 154 (1935). 
H. Feshbach, Ann. Phys.(N.Y.) 5, 357 (1958). 
H. Feshbach, Ann. Phys.(N.Y.) 19, 287 (1962). 
E. Tiesinga, B. J. Verhaar, and H. T. C. Stoof, Phys. [15] 
Rev. A 47, 4114 (1993). 

C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. 
Mod. Phys. 82, 1225 (2010). 

J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. [16] 
Rev. 106, 162 (1957). 

C. Regal, C. A.and Ticknor, J. L. Bohn, and D. S. Jin, [17] 
Nature 424, 47 (2003). 

P. Wolfle, Reports on Progress in Physics 42, 269 (1979). [18] 
T.-L. Ho and R. B. Diener, Phys. Rev. Lett. 94, 090402 
(2005). 

C.-H. Cheng and S.-K. Yip, Phys. Rev. Lett. 95, 070404 [19] 
(2005). 

V. Gurarie and L. Radzihovsky, Ann. Phys. 322, 2 [20] 
(2007). 

C. A. Regal, C. Ticknor, J. L. Bohn, and D. S. Jin, Phys. 



Rev. Lett. 90, 053201 (2003). 

J. Zhang, E. G. M. van Kempen, T. Bourdel, 
L. Khaykovich, J. Cubizolles, F. Chevy, M. Teichmann, 
L. Tarruell, S. J. J. M. F. Kokkelmans, and C. Salomon, 
Phys. Rev. A 70, 030702 (2004). 

C. H. Schunck, M. W. Zwierlein, C. A. Stan, S. M. F. 
Raupach, W. Ketterle, A. Simoni, E. Tiesinga, C. J. 
Williams, and P. S. Julienne, Phys. Rev. A 71, 045601 
(2005). 

C. Ticknor, C. A. Regal, D. S. Jin, and J. L. Bohn, Phys. 
Rev. A 69, 042712 (2004). 

J. P. Gaebler, J. T. Stewart, J. L. Bohn, and D. S. Jin, 
Phys. Rev. Lett. 98, 200403 (2007). 

J. Fuchs, C. Ticknor, P. Dyke, G. Veeravalli, E. Kuhnle, 
W. Rowlands, P. Hannaford, and C. J. Vale, Phys. Rev. 
A. 77, 053616 (2008). 

S. T. Thompson, E. Hodby, and C. E. Wieman, Phys. 
Rev. Lett. 95, 190404 (2005). 

Y. Inada, M. Horikoshi, S. Nakajima, M. Kuwata- 
Gonokami, M. Ueda, and T. Mukaiyama, Phys. Rev. 
Lett. 101, 100401 (2008). 



11 



R. A. W. Maier, C. Marzok, C. Zimmermann, and P. W. [39 
Courteille, Pliys. Rev. A 81, 064701 (2010). 
M. H. Szymaiiska, B. D. Simons, and K. Burnett, Phys. [40 
Rev. Lett. 94, 170402 (2005). 

M. H. Szymaiiska, K. Goral, T. Kohler, and K. Burnett, [41 
Phys. Rev. A 72, 013610 (2005). 

J. R. Taylor, Scattering Theory (Dover, New York, 2006). 

B. H. Bransden and C. J. Joachain, Physics of Atoms [42 

and Molecules (Prentice Hall, 2003). 

R. G. Newton, Scattering Theory of Wave and Particles 

(Dover, New York, 2002). [43 

A. J. Moerdijk, B. J. Verhaar, and A. Axelsson, Phys. 
Rev. A 51, 4852 (1995). [44 
K. B. Gubbels and H. T. C. Stoof, Phys. Rev. Lett. 99, 
190406 (2007). [45 

B. Gao, J. Phys. B: Atomic, Molecular and Optical 
Physics 37, 4273 (2004). [46 
F. Chevy, E. G. M. van Kempen, T. Bourdel, J. Zhang, 

L. Khaykovich, M. Teichmann, L. Tarruell, S. J. J. M. F. [47 
Kokkelmans, and C. Salomon, Phys. Rev. A. 71, 062710 
(2005). [48 

D. E. Manolopoulos, J. Chem. Phys. 85, 6425 (1986). 

H. T. C. Stoof, J. M. V. A. Koelman, and B. J. Verhaar, [49 
Phys. Rev. B 38, 4688 (1988). 
S. Falke, H. Knockel, J. Friebe, M. Riedmann, E. Tie- [50 
mann, and C. Lisdat, Phys. Rev. A 78, 012503 (2008). [51 
A. L. Fetter and J. D. Walecka, Quantum Theory of 
Many-Particle Systems (Dover, 2003). [52 
L. Austen, PhD thesis (University College London, Lon- 
don, 2011). [53 
M. Holland, S. J. J. M. F. Kokkelmans, M. L. Chiofalo, 
and R. Walser, Phys. Rev. Lett. 87, 120406 (2001). [54 

E. Timmermans, K. Furuya, P. W. Milonni, and A. K. 
Kerman, Phys. Lett. A. 285, 228 (2001). [55 
J. Levinsen, N. R. Cooper, and V. Gurarie, Phys. Rev. 
Lett. 99, 210402 (2007). 



V. Gurarie, L. Radzihovsky, and A. V. Andreev, Phys. 
Rev. Lett. 94, 230403 (2005). 

K. E. Strecker, G. B. Partridge, and R. G. Hulet, Phys. 
Rev. Lett. 91, 080406 (2003). 

J. CubizoUes, T. Bourdel, S. J. J. M. F. Kokkelmans, 
G. V. Shlyapnikov, and C. Salomon, Phys. Rev. Lett. 

91, 240401 (2003). 

S. Jochim, M. Bartenstein, A. Altmeyer, G. Hendl, 
C. Chin, J. H. Denschlag, and R. Grimm, Phys. Rev. 
Lett. 91, 240402 (2003). 

R. A. Barankov and L. S. Levitov, Phys. Rev. Lett. 93, 
130403 (2004). 

R. A. Barankov, L. S. Levitov, and B. Z. Spivak, Phys. 
Rev. Lett. 93, 160401 (2004). 

A. V. Andreev, V. Gurarie, and L. Radzihovsky, Phys. 
Rev. Lett. 93, 130402 (2004). 

E. Altman and A. Vishwanath, Phys. Rev. Lett. 95, 
110404 (2005). 

B. Dobrescu and V. Pokrovsky, Phys. Lett. A. 350, 154 
(2006). 

L Tikhonenkov, E. Pazy, Y. B. Band, M. Fleischhauer, 

and A. Vardi, Phys. Rev. A 73, 043605 (2006). 

S. Matyjaskiewicz, M. H. Szymahska, and K. Goral, 

Phys. Rev. Lett. 101, 150410 (2008). 

V. Gurarie, Phys. Rev. A 80, 023626 (2009). 

M. Iskin and C. J. WiUiams, Phys. Rev. A 77, 041607 

(2008). 

C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 

92, 040403 (2004). 

M. Jona-Lasinio, L. Pricoupenko, and Y. Castin, Phys. 
Rev. A 77, 043611 (2008). 

J. Levinsen, N. R. Cooper, and V. Gurarie, Phys. Rev. 
A 78, 063616 (2008). 

C. Lovelace, Phys. Rev. 135, B1225 (1964). 



